X 



Some peculiarities of transition from discrete 
to continuum Fourier series in lattice theories 



(N 
O 

o . 

CN ■ Vladimir K. Petrov 

O 

CD 

Q 



A^. A^. Bogolyubov Institute for Theoretical Physics 

National Academy of Sciences of Ukraine 
252143 Kiev, Ukraine. 11.12.2002 



> 
00 

Q ■ Abstract 

(N ; 

Transition from discrete to continuous Fourier series is studied for 

CN ■ the functions becoming singular in the transition. Conditions are spec- 

. ified when summing replacement by integration is inadmissible. 

> . 

Oh! 1 Introduction 

D 

Fourier series in discrete variables, or discrete Fourier series, is a routine 
instrument of computation on finite lattices with cyclic border conditions. 
Quite often it becomes necessary to remove lattice regularization with lat- 
5^ \ tice length N approaching infinity. Ordinarily lattice theory operates with 

regular functions. However, these functions may include terms which, in 
fact, play the role of a regularization parameter. In other words, even when 
such functions are analytical in the whole physical region for any finite N , 
the transition to — > 00 may convert regular functions into distributions. 
For the function with such nonstandard behavior, ordinary replacement of 
summing by integration for N ^ 00 may appear inadmissible. 

In this paper we focus our attention on a frequently occurring case, when 
such "regularization" terms are inversely proportional to the lattice length. 
The conditions are specified when such replacement may or may not be done 
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and instances are discussed, when this replacement is definitely erroneous 
and relevant corrections are needed. We also suggest two different methods 
for correction computation. 



2 Some Remarks on Discrete Fourier Series 

In the previous work [|l| we considered the consequences of cyclization pro- 
cedure 

oo 

F{^)= J2 ^ ('^ + 2™) (1) 

n=— oo 

application to the function r {^p) of the continuous variable (p. Similar proce- 
dure may be applied to the function Fk of a discrete variable k = 0, ±1, ±2, .... 
We assume that is a tempered function, i.e. that there exists a constant 
c < cxD so that 

\Fk\<\k\''. (2) 

for \k\ ^ oo 

To simplify the problem, we associate Fk with some periodic distribution 

F{^)=F+{^ + ie)-F_{if-te) (3) 

where e is a routine positive infinitesimal parameter and F± (ip) are analytical 
functions in upper/lower complex half-plane 

oo —1 

F+ (^) = ^ Fk exp {ik^} ; F_ ((^) = - exp {ik^} (4) 

k=0 k=—oo 

Bearing in mind that a distribution is tempered, if and only if it is a finite 
order derivative of some continuous tempered function and that any tempered 
distribution may be presented in the form (j^) (see e.g. 0, ^, ^). Moreover, 
every periodic distribution is tempered 

As it is seen from (D the expression coincides with an Abel-Poisson 
regularization of Fourier series 



oo 



F{^)=reg Fkc''"' = Fkc''''-'^'^ (5) 



k=—oo k=—oo 
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Since Fourier coefficients Fj. in (j^) coincide with the considered discrete 
function, they are related to F by the standard equation 

^k = Tr\ F (ip) exp {-iifik} dip (6) 

where an arbitrary real number a is usually chosen as a = or a = tt for 
convenience reasons. 

Now one can easily check that applying cyclization procedure to F^ in 
with the period 

oo 

Fk = ^ Fk+vN- (7) 

11=— oo 

and making allowance of the Poisson relation 

oo oo 

J2 e""^ = 27r J2 2vrr) (8) 

k=0 r=—oo 

one gets 

= ^ £ j ^^F{^)5[lp~ exp {-i^k} dip (9) 

From d^) we see that a sufficient condition for Fk to exist is the existence 
of the product of distributions F [ip] and 6 {ipN — 27m). Unfortunately, for 
an arbitrary distribution (in our case for the arbitrary F (ip)) such product 
is not defined ^, ^, which requires a more detailed consideration. 

It is known that Dirac 5-function is rigorously defined as functional {6 {p>) , {p>)) = 
(/) (0) on the space C° of continuous functions 4> {p>) G C° @]. Therefore, if 
F {p>) G then F {p>) (f) [ip) G C", i.e. F (ip) is a multiplicator in C° so for 
such F {(p) the considered product is well-defined. 

It is hardly worth mentioning that since F [ip) is given by the Fourier 
series (^) this function is not absolutely arbitrary and had to obey specific 
conditions. However, even regarded as too restrictive, Dirichlet conditions 
require function F {(p) only to be piecewise monotonic, piecewise differen- 
tiable and bounded for < (y9 < 27r (see e.g. 0). Under given conditions 
Fourier series (|]) converge to F (ip) at all points of [0, 2tt] , but at some points 
p>i, of this segment it may have first order discontinuity. Despite the fact that 
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F [ip) is defined in all these points as F (ip^) = \ (^F [ip^ — 0) + F [ip^ + 0)^ 

it is no remedy for the discontinuity. Therefore, when such discontinuity 
points coincide with ip = 2TTn/N, special consideration is necessary. The 
distribution theory allows to consider even more singular function F [ip). 

If F {(p) is singular, then function Fk defined in @ may appear to be 
nonperiodic, i.e. F^ 7^ ^^/c+at. This case is very similar to one considered in 
^ where the procedure (P was applied to a function, the Fourier integral of 
which has singularity. Therefore we confine ourselves to a simple case and 
take an example when F {(p) has the simple pole 

Fi^) = . (10) 

1 — exp {lip} 

From we get 

= ^ £ l-exp{.y.} ^ - ^) ^-^^^^ ^''^ 



Without loss of generality we may take as an odd number 

N J^^ -1 + ^ip> + ^iip>^ + O (v3^) ip 



JV-1 



1 ^ / expl-^gf^A:} expj^g^fc} 
iV^ I l-exp{-z^} l-exp{i^} 



n=l 



Making allowance of the fact that for any function / {ip>) regular and nonva- 
nishing in the vicinity of (/? = we may write 

/ -5 (^) = -/ (^) 5' = r (0) 5 (y.) , (13) 



we easily get from (|T2| 



AT-l 



Fi 



f->7 , -1 1 2 2iTnk 27rn(l— fc) 

2k + l ly. COS COS 



k 



N iV ^ 2 sin 

n=l J» 



It is evident that the second term in the right part of the equality ([I^ 
is invariant under substitution k ^ k + N . Such term appears as the sum 
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of contributions of the points = 21171/ N] n ^ the first term being the 
contribution of the pole of function F {ip) at = (i.e. n = 0). This term 
is evidently nonperiodic. 

When function F {(f) is bounded at all ip = 2iTn/N and continuous in the 
vicinities of those points, we regard that the sufficient condition for F^ to be 
periodic is fulfilled. Indeed, from (|) we obtain 

1 ~ f27in\ ( 2nn,] 

(15) 

Therefore, since Fk is presented as the finite Fourier series it is an evidently 
periodic function with the period N. It is clear that F (^^), as a function 
of discrete variable n, is periodic with the same period by definition. 
Making allowance of the simple identity 



1 r 2nk 

— > exp < —i——n 

k=0 



E ^n"" (16) 



one may write the transform inverse to (p!5|) 

J2F,exp\-z—-k\ (17) 



k=0 



We see that ([T7|) is nothing but the discrete version of Fourier series 



as well as (|T5D is the discrete version of (^. 

It is clear that the above condition on F {(f) may be related to the con- 
vergence of series (|]) or rather to the behavior F^ for k — ±00 that, in its 
turn, is related to the convergence of series in (^). For instance, necessary 
conditions on F {(p) are fulfilled when series (|]) is uniformly convergent on a 
closed segment — vr < (f> < n (see e.g. [0). 

However, one may impose on Fk a weaker restriction and consider F {(p) 
which has a limited variation, i.e. such functions which for arbitrary partition 
— TT < (fy < (f>v+i < TT and arbitrary obey 



N-l 



v=0 

for some M < cxo. 



E F(v9,)-F(v9.+i) <M. (18) 
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A function with a limited variation in a segment is obviously bounded 
on it. Moreover, as it is shown by N.Weiner (see e.g. such function is 
continuous, if (and only if) 



This condition, however, may be regarded as too restrictive as well, just 
for that very reason that it claims 'good' behavior of F {ip) for all cp, whereas 
it is actually needed only at infinitesimal vicinities of = 2TTn/N. Hence 
there is no reason to expect that each particular F {p) that is presented as 
discrete Fourier series ([T7|) obligatory obeys such conditions . 

3 Regular izat ion Removal and Transition to 
Continuum 

The discrete Fourier series (|^) and (|1^) are commonly used in lattice theories 
with periodic boundary conditions. Quite often one needs to apply a proce- 
dure inverse to (|^, which involves the removal of the lattice regularization 
and the lattice length then tends to infinity. 

This section considers peculiarities of the transition from discrete to con- 
tinuous Fourier series with N ^ oo. As it is said above, we are interested 
in a very particular aspect of such a fundamental problem as transition to 
continuum in a lattice gauge theory, and namely in the transition to contin- 
uum of the discrete Fourier series for functions which are regular at finite 
A^, but convert into distributions with N ^ oo. Naturally, some impor- 
tant issues of transition to continuum are outside the limits of this study, 
nonetheless, the suggested approach allows, in particular, to compute the 
fermion contribution on an extremely anisotropic lattice §]. 

Let F (ip) in (|I5]) obey the conditions of Lebesgue theorem, i.e. F (p) is a 
bounded function with a set of measure zero discontinuity points. In this case 
F {(fi) is Riemann integrable, i.e. the difference between upper and lower sum 
Darboux vanishes for arbitrary partition (— vr < ipv < ^v+i < tt; u = 0,1, N — 1) 
when the number of terms A^ infinitely increases^. For N oo this allows 

"'^For piecewise-monotonic F (ip) the sum in (|l5| ) decomposees into partial sums over 
monotonia pieces. Each partial sum coinsides either with upper or with lower sum Dar- 
boux. 




(19) 



k=- 



K 
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us to make a change in ( |T3| ) 



27rn _ 1 1 
N 27r 

n=0 



1 ' , 1 /"^'^ 



which transforms (|15D into (|]). 

For any finite value s we may put s/iV — >■ 0. If, however, s = is a 
singular point of function F {(p) we need to preserve the information about 
the direction in which s/N approach zero in a complex s-plane 

s/A^ ^ e exp {i arg s} (21) 

When, however, a transition for Fk in (|l^ is of interest, one needs Fk to 
obey the conditions of Lebesgue theorem. If it is the case, then transition to 
continuous variable 2iTk/N ip may be fulfilled by the following substitution 

1 1 C'^'^ ~ ~ ~ /27r \ 

To estimate the error, introduced by the procedurefpOf) we apply the 
Poisson formula to truncated functions 

F(y,) = p(^)' 0<^<27r; ^23) 
1 0; otherwise. ^ ' 



that gives 

= E ^ / F{^)e.^{-i^kN}d^ (24) 



n=0 ^ ' k=—oo 

SO, taking into account the definition (|^) we may write 



n=0 ^ ^ fc=-oo 



It is clear that the application of procedure ( ]20| ) to the left part in ( pS]) 



gives but the term Fq. Therefore, to obtain an exact result, one had to 
include the correction A which is defined as 

oo 

A = ^ FfcAT = ^ {FkN + F_kN) (26) 

k=iO k=l 
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For A to vanish, it is sufficient that for all k 



F exp {-iipkN} ^ 0; N ^ oo 
Since kN is an integer number we may rewrite (^) as 



(27) 



F (ip) exp {-i - ao) kN} 0; N ^ oo 



(28) 



where (Tq = 27mo/N and no is an arbitrary integer number. In other words, 
the result of application (^) to left-hand member of (0) coincides with 
right-hand one only if the condition (^Sf) is satisfied for any cxo and all k 0. 

According to Riemann-Lebesgue lemma (see. e.g. 0]) all F^n in (^) with 
k ^ turn into zero with — oo, if either the integral in (p4|) converge 
absolutely for k = or F (</?) has limited variation on < (/? < 27r. 

For singular enough F{(p), the conditions of Lebesgue theorem may be 
violated and procedure ( pOD transforms the sum into improper integral, which 
value will depend on the choice of the way singularity is approached. The 
result of proceeding to limit in expressions like (^) was studied in detail 
IPHI , Various type singularities of F{ip), located infinitely near to ip = 
were considered. Results obtained in [|10|, |lT| allow to conclude that the 
condition ( PSj) is satisfied, if F [ip) is regular at infinitesimal vicinities of all 
if = 27m/ N, where n = 0, — 1. 

Commonly the function F {(p) in ( PS] ) doesn't depend on A^. In a lattice 
gauge theory, however, F {(f) may include small parameters, such as fermion 
mass term ma, which tends to zero with A^ — oo, provided that inverse 
temperature 1/T = aN remains finite. Actually, it introduces a dependence 
on A^ into F {ip). Therefore, it is not without interest to consider a case when 
F (ip) includes some 'regularization' term which is proportional to 1/A^, i.e. 
singularity of F (ip) is located in a complex ip plane at finite distance ips ~ 
i/N from the real axis. In this case the increase of exponent oscillations in 
(|28| ) is matches singular point ips approaching physical region. In Appendix 
we suggest a simple recipe to modify the results obtained in []TU|, for this 
case. 
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4 Examples of nonstandard behavior 

Let us consider an example of a periodic function with the simple pole 
27m \ (.27m xl\ ^ f Z' .2nn x 



= l-exp -- exp -™ U— --)^ (29) 



It is clear that the role of regularization term in (p9D is played by x/A^, where 
X is finite and Im x = 0. 
Writing series expansion 



oo 



= ^(-x) +^ f^(x)e'^^^ -^(-x)e-'^^^j (30) 



r=0 



and collecting coefficients at r = modulo A^, we obtain 

27rin-x\~l „ , ^ 9 { K) X — ^ i,27rin-x 9{ — x) X — ^ _^27rin- 



/ ^ ^ 1 - ^ 1 - ^ 

/c=0 A;=0 

(31) 

As it follows from ([TBI) one may write 

n=0 



for nonnegative integer m < N and n < N. Therefore, from (^Tj) we finally 

get 



1 -1 1 

-L V ^ / 27Tin — >c \ 27vin — >c J. , ^ , „™ ^ 

]vE(l-^^) e-^ = ^-^-^^(-x)(l-5o"^) (33) 



n=0 



Proceeding to limit A^ ^ oo in ( p3| ) is trivial, because the right-hand member 
doesn't depend on A^. 

On the other hand, the application of procedure ( pO]) to the left-hand 
member in (^3)) gives 



exp{ zv^m} -rfy, = signum(x) + 0(-x)5^ (34) 



27r 1 — exp {i(p — e signum (x)} 
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We see that the procedure leads to a result which differs essentially 
from an exact one given in (pSl). The reason is that the difference between 
upper and lower sum Darboux computed for F{2mn/N) in (^) doesn't 
vanish for N ^ oo, thus the sum in (|3^) doesn't convert into Riemannian 
integral when N ^ oo. It may be anticipated, because Riemannian integral 
is defined only on bounded function, whereas (pSj) remain finite only for finite 
N. 

Let us compute now the correction A defined in (^) . In considered case 
FkN is given by 

'•^^ exp{-i^kN} dip 







in 



l-exp{z(^-f }27 

= [e (x) e{k) + e (-x) - e (-x) e {-k)] exp {-kx} (35) 

and we easily get 

A = e (x) (36) 



Thereby it is confirmed that the correction (|36D added to (|34D allows to obtain 
the exact result (RBf). 



Applying the differentiation with respect to a parameter x one may easily 
repeat the computation for arbitrary rational function of exp {— 27rm/A^}. 
Previously considered case may be easily generalized in another way, if we 
consider 

r°° 1 

F{^)= f (x) J-. ^rfx (37) 

J-00 l-expl^v^-^l 

where the integral has to be interpreted in the sense of a Cauchy principal 
value. Making allowance for (pSf) we get 



(38) 

On the other hand, the procedure ( pO]) applied to the left-hand member 
in (|38D leads to 

/(x)rfx \ .^rnd^ 



J 1 — exp {i(p — e signum (x) } / 27r 

POO pO 

/ /(x)rfx-(l-O / /(x)rfx (39) 

^0 J-00 
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which disagrees, as it may be anticipated, with 

Let us now consider as an example the periodic function with logarithmic 
singularity. From ( [|13|6. 1.1(10)) one can easily get 



N-l , r r, ■ 



n^l-expl^jj^l-e- (40) 

that leads to 

X;i„(l-exp|?^|)=ln(l-e-) (41) 

n=0 ^ ^ 

On the other hand, the procedure (^) gives 

1 r 

— / In (1 — exp {i {if + ie signum (x))}) dip = Q (42) 
27r 



that evidently disagree with (^Tj) at finite x. 

In some instances the procedure (|^) leads to ambiguous results and be- 
comes inadmissible. Indeed, from (|13[4.4.5(7)) one can get 

(.) ^ 1 f (^'^'") ) e--^ ^ ''"-/^) . (43) 

Hence, one obtains a well defined result for N ^ oo and any finite m 

lim Qn (x) = g (x) = ^— . (44) 

N^oo COS X 

On the contrary, the procedure (PD|) gives 



1 r e*^^cos(|m)sin(^) , , 

Q (x) = lim - / TT^^S^^^ 

A^-^oo 27r cos (rj) — cos (-^ j 

where rj = e exp {i arg x}. 

Computation in (^) of Q (x) leads to an ambiguous result, because, on 
the one hand, exp {iipN} turns integrand into zero for all ip except the vicin- 
ity of = 0, on the other, - the integrand has the pole aX (p = 27] that 
is located infinitely near to = 0. Therefore, Q (x) includes the factor 
lim^v^oo exp {2i?7A^} which cannot be defined unambiguously without addi- 
tional assumptions. 
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5 Appendix. Asymptotic Series for Distribu- 
tions 

Here we briefly sketch some of the results obtained in [1^ for the asymptotic 
expansions of exp {'^iipN} F (ip) for N ^ oo with necessary modifications 
which we suggest for a case when the regularization term in F {ip) is propor- 
tional to 1/A^. 

As it is known (see e.g. 0), temperate distribution is defined as func- 
tional on infinitely differentiable fast decreasing test functions $ {(p) i.e. 
lim|^l^oo Iv^P^*-^-* (p) = 0; for any m > 0; p > 0. As it pointed out in |TD[, 
one may write the asymptotic (Poincare) series for the functional defined on 
any test function $ {p) 

T^^ I exp {TipN} F {p) $ {p) dp c^Y.^^ {N) I (^) $ (^) dp^ (46) 

n=0 



where ijj^ (x) for x ±oo constitute asymptotic sequences obeying 



So, differentiating n-times (|46|), we may get in particular 

p^C^ {p) = (48) 

that means 

n-1 

iv) = E C'nV^'^ iv) (49) 



fc=0 



The simplest form of C„ (p) ip^ (N) satisfying both(^9|) and condition (^Tf) 
is 



n— 1 n— 1 

{V) (N) = J2 C^ki^n (N) ip) ~ Pnk (In N) N-'^-H^'^^ (p) 

k=0 k=0 

(50) 

where polynomials (In A^) and constant r are defined by the specific form 
ofF{p). 
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Bearing in mind that any tempered (not obligatory periodic) distribution 
may be presented in the form analogous to (Q) (see e.g. [Q) where F± (y?) 
are analytical functions in upper/lower complex half-plane if 

F{^)=F+{^ + ie)-F_{^-ie) (51) 



it can be shown that ||T0| 

F± ± ie) exp {±i^N} ~ 0, ^ oo (52) 

and 



F± {if ± ie) exp {Ti^N} ^ ^ (A^) ^^"^ (</^) , N oo (53) 

n=0 

As already mentioned, we are interested in a specific case, when singular- 
ity of F± is located not at infinitesimal distance iie but at finite^ one -^ix/N 
from the real axes of ip and approaches the real axes only with N ^ oo. 
Thereby, instead of a definition given in (^) we may write 

T^^ I exp {Tt^N} F± ± z^) $ (^) (54) 

and after a simple shift (f ^ (f^i {x/N — e) we obtain 

(55) 

The test functions $ {ip) and their derivatives {ip) should not be oblig- 
atory analytical, but they are infinitely differentiable. Hence, their Taylor 
series may be taken as asymptotic one, no matter it is convergent or divergent 
in the area in question [|1^ 

m=0 ^ ^ 

(56) 

■^In this section we assume >c > Q. 
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Taking into account that parameter e is infinitely small in comparison 
with x/N we obtain from (^) 



^ — n ^ — n \ / \ 



71=0 m=0 

and for — oo we get 



n=0 

or 



n=0 



that with (|53|) and ( ^21) gives after change — > /cA^; x ^ kx 

( k 

exp {=Fiv5/cA^} F± ± i — 



exp{-fcx}Er=o^n /or fc>0 

for k <0 



We see that correction A in ( pS] ) may be written in a form 

A = 5^ exp {-xk} — 5^ C± (fcAT) 5W _ ^o) ciy. 

fc=l '^'J n=0 



27r 

fc=i 



(57) 



oo 

T^i^ = e-^Y.^^ (^) (0) (l + O (^)) (58) 



exp {Tz^AT} F± ± ^ exp {-x} (A^) ^^"^ (^) (59) 



(60) 



^ oo 

-Y^exp{-xk}C^{kN) (61) 



As it is seen from ( |6TD the correction A is presented as a series which 
is a discrete version of Laplace transform, called Z - transformation (see 
e.g. [T^). Since we consider x/N — » +0, coefficients (x) in this series 



are the same as those computed for the asymptotic expansions in |12]. It 



should be noted, that distributions studied in |]T2[ are nonperiodic. Hence, if 



direct application of the method developed in ||I2[ to the periodic distribution 
under consideration is technically difficult, it may be convenient to apply 
the cychzation procedure considered in |]1|| to those distributions that have 
already been studied in [Q. 
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6 Conclusions 



Some particularities of the transition from discrete to continuous Fourier 
series in a limit N ^ oo are discussed. A case is studied when the regular- 
ization parameter is proportional to and the considered functions are 
analytical for N < oo, but they are converted into generalized functions or 
distributions when regularization is removed, i.e. N — >■ oo. 

Some simple and rather general examples of functions with 'nonstandard' 
behavior were presented to clarify the conditions under which ordinary re- 
placement of summing by integration is inadmissible and the corresponding 
corrections A were computed. It is also pointed out that when the peri- 
odic function F [if) obeys the conditions of Riemann-Lebesgue lemma, the 
corresponding correction turns into zero. 

We suggested a modification of asymptotic series expansions proposed in 
|llO| , |Tl| for a case with regularization parameter proportional to and 
applied the modified series for the computation of the correction A. 
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